Are defect models consistent with the entropy and specific heat of glass-formers? 
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We show that point-Uke defect model of glasses cannot explain thermodynamic properties of glass- 
formers, as for example the excess specific heat close to the glass transition, contrary to the claim 
of J.P. Garrahan, D. Chandler [Proc. Natl. Acad. Sci. 100, 9710 (2003)]. More general models 
and approaches in terms of extended defects are also discussed. 



Kinetically constrained models (kcm) display trivial 
thermodynamics but non trivial dynamics Q, y, y| . They 
encapsulate in a specific way the 'free volume' idea ac- 
cording to which the slow dynamics of glasses is due to 
rare, localized mobility defects 0, 01 . For most kcms, 
slow dynamics is indeed due to defect motions 3]. A 
remarkable aspect of these models is that although the 
dynamical rules are very simple and local, the emergent 
defect dynamics lead to highly non trivial physical be- 
haviour. As a consequence kcms provide a theoretical 
framework for understanding many of the remarkable dy- 
namical properties of fragile glass-formers, such as super- 
Arrhenius relaxation time i^, iTj, dynamical hetero- 
geneities H, 0, 0, ^3 ^'^d viscosity/diffusion decoupling 
|l2]. That these models may in fact describe quantita- 
tively the physical properties of glasses has been strongly 
advocated in a series of papers by J.P. Garrahan and D. 
Chandler (gc) and investigated further in jl^]. 

One of the major tenet of kcms is the complete decou- 
pling between dynamics and thermodynamics, at vari- 
ance with the traditional view of Adam and Gibbs , 
and others 0, 0| , in which the decrease of configura- 
tional entropy with temperature is the fundamental un- 
derlying mechanism for the viscous slowing down. It is 
therefore a priori surprising that kcms could have any- 
thing sensible to say about the entropy and specific heat 
of real materials - even if these turn out to be accept- 
able models of their dynamics. Nevertheless, the claim 
made in is that these models are also in quantitative 
agreement with the specific heat of these materials. More 
precisely, assuming a perfect gas of free defects, GC give 
a formula for the specific heat at fixed pressure per par- 
ticle of the liquid (in excess of that of the solid), which 
we reproduce from ref . |lO| to be: 



where J\f is the number of physical molecules contributing 
to enthalpy fluctuations in a coarse-graining cell, J is the 
enthalpy cost, in kelvins, for creating a mobile cell (re- 
member that liquids are considered here at constant pres- 
sure, hence the importance of distinguishing between en- 
thalpy and energy) . The concentration of mobile cells c is 
assumed to be given by c(T)/(l — c(T)) = ylexp(— J/T), 
where the numerical prefactor A is argued by GC to be 
rather large: In A is akin to an entropy gain Asq/^b asso- 
ciated with the creation of a mobile cell. Note that there 
is a small difference between the specific heat in excess 
to that of the glass (considered in 10]) and that in excess 
of the crystal; this difference is however irrelevant for the 
present discussion. 

In the following we first adress the results of show- 
ing that even if the value of ACp at Tg seems to match 
experimental values, the temperature dependence of ACp 
predicted by Eq. is in total disagreement with exper- 
imental values. We then discuss in detail the derivation 
of ACp(r) and the physical assumptions behind Eq. 
This derivation leads to a formula different from Eq. 
and makes clear that the defect contribution to ACp(r) 
is in fact completely irrelevant close to the glass tran- 
sition. In order to avoid a contradiction, one has to as- 
sume that the main contribution to both the entropy and 
the specific heat comes from the immobile regions. Al- 
though this is most reasonable physically, it implies that 
the models considered by GC are incapable of reproduc- 
ing the thermodynamics of fragile glasses (and actually 
neither of strong glasses). 

A clean way to test Eq. Q , that gets rid of any ambi- 
guity in prefactors (we in fact disagree with the prefactor 
JV in Eq. see below), is to rewrite it as: 
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This formula is expected to be valid for small c 
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Aexp(— J/T), where the description of glassy dynamics 
in terms of rare, dilute, defects could make sense - cor- 
responding to a large enthalpy of creation for mobility 
defects, J/Tg > 1. 

Using viscosity data at fixed (atmospheric) pressure^ 
GC are able to estimate the value of J/Tg, found to 
be equal to 16.7 for 3-bromopentane (3BP) and 26.7 
for ortho-terphenyl (OTP), indeed quite large compared 
to unity. Fig. 1, on the other hand, gives a plot of 
ACp{T) / ACp{Tg) , obtained from Eq. ISJ, as a function 
of T/Tg for J/Tg = 10, 20, 30 and compares it to the ex- 
perimental values for 3BP [13 and OTP Jl^]. (Note the 
log scale on the y-axis.) It is clear that Eq. (jSJ is totally 
incompatible with the data when J/Tg is large: it varies 
much too fast with temperature, and increases with T 
instead of decreasing (a similar criticism has been made 
by C. Angell and coworkers for excitation models). 

An acceptable fit of the data with Eq. ^ requires 
J/Tg « 1.3, meaning that within this description the 
density of defects is not small at all - but this of course 
is then incompatible with viscosity data. For consistency, 
Eq. ^ cannot be used for temperatures such that c ~ 1 ; 
GC restrict their analysis to a temperature range between 
Tg and T1/2, where T1/2 is defined as the temperature at 
which the logarithm of the viscosity expressed in Poise, 
or of the relaxation time expressed in picoseconds is half 
its value at Tg (see Fig. 1). One can clearly see that, 
even in this restricted interval of temperature close to Tg 
(for which c is of the order or less than 10^^ ^^'^ 
theoretical prediction is in complete disagreement with 
the experimental data. 

The previous discussion suggests that the thermody- 
namic predictions of the models considered in are in- 
compatible with experiments. As we discuss below, there 
is in fact no real contradiction: the order of magnitude of 
the defect contribution to the specific heat is simply much 
smaller than the experimental value between Tg and Ti^2: 
which is instead dominated by the thermodynamic prop- 
erties of immobile regions. In the models considered in 

llfl||. however, the immobile regions give a trivial con- 
tribution to thermodynamics properties. In order to be 
consistent, one should assume that these regions domi- 
nate the thermodynamics but are irrelevant dynamically. 

Let us now reconsider in detail the derivation of Eq. 
taking explicitely into account that experimentally both 
the specific heat and the excess entropy are found to be of 
order ks per particle or larger at the glass temperature - 
indeed a very important "stylized fact" of fragile glasses. 
An immediate consequence of this experimental fact is 
that the excess entropy of the liquid cannot be due to 
point-like defects only: since these must be dilute to lead 
to large relaxation times, their contribution to the excess 
entropy must be small. Within a simple coarse-grained 
defect description, the total excess entropy per cell (with 
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FIG. 1: Normalized excess specific heat as given by GC's for- 
mula, Eqs. Q and (|5J, for J/Tg = 10,20,30, compared with 
the experimental data for OTP and 3-BP, for which GC esti- 
mate, from viscosity data, J/Tg = 26.7 and 16.7, respectively. 
Note the log scale on the j/-axis. The two vertical lines corre- 
spond to Ti/2 defined as the temperature at which the loga- 
rithm of the viscosity expressed in Poise, or of the relaxation 
time expressed in picoseconds is half its value at Tg. 



respect to the crystal). As, is: 

As = -kB[cAnc+ (1 - c)ln(l - c)] + cAsq + As^, (3) 

where the first term corresponds to the entropy of the 
ideal gas of defects, the second term to the entropy dif- 
ference between mobile and immobile cells (see above), 
and the last term is the configurational entropy asso- 
ciated with frozen cells In order to get the ex- 
perimental order of magnitude of As at Tg and at the 
same time have ACp dominated by defect properties, one 
must assume that (i) Asc is of order J\fkB, hence much 
larger than the defect entropy itself; (ii) Asc must then 
be temperature independent, otherwise it would give the 
leading contribution to ACp. Although large, the ex- 
cess entropy Asc should not play any physical role in the 
slowing down of the material in this pictu re ( in strong 
contrast with the view expressed in [Mlialia)- There- 
fore, the excess specific heat at constant pressure per cell., 
Acp(T) = TdAs{T)/dT, is given by: 



Acp{T)^kB-T^^ 



kBiJ/Tfc{T), (4) 



where we have used c{T) « exp(Aso — J/T), valid in 
the limit c <g; 1. Note that the functional form of c(T) 
cannot be derived within the model and it is in principle 
obtained from the coarse-grained procedure that should 
map glass-forming liquids to kcms. Since until now no 
such procedure has been developed one has to postulate 
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an a priori functional form c(T). The simple physical 
assumptions described above correspond to the one as- 
sumed in [loj . However, we will show below that our 
conclusion is indeed independent of the form of c(T). We 
have obtained a specific heat per cell given by ksiJ /T^c, 
and therefore a specific heat per molecule given by: 

^C,{T)^kB[pj\{T)N"\ (5) 

a factor smaller than Eq. . The same expression is 
obtained by looking directly at the enthalpy fluctuations 
|22l | . Taking very favorable values for the parameters c = 
0.01, J/Tg = 30 and A/" = 10 leads to ACp(Tg) « O.lfcs, 
already a factor 100 too small compared to experimental 
values (9fcB for 3-BP, 13fci3 for OTP). The only way to 
avoid a contradiction is that the temperature dependence 
of Asc is strong, and that in fact all the thermodynamics 
is contained in that very contribution. 

The above conclusion is based on comparing the exper- 
imental specific heat to that predicted by defects only. 
A more general argument (see also 0, H^l), indepen- 
dent of the shape of c(T), relies on the entropy change 
between - say - Tg and T1/2, which is experimentally 
found to found to be of order fc^ (« Aks for OTP). 
Thus, not only one finds that Asc(r) has to be much 
larger than (— fc_B[clnc -I- (1 — c)ln(l — c)] -I- cAso) at 
Tg, as discussed previously, but also that Asc(T) has 
to provide the leading contribution between Tg an d Ti 
since in this regime the ideal gas entropy of defects Q is 
negligible As a consequence, using the relationship 
Acp(T) = TdAs{T)/dT, one finds again that the main 
contribution to ACp comes from the immobile regions 
close to Tg. Note that this result is independent of the 
functional form c(T), as long as defects are dilute, c ^ 1, 
a necessary condition for the GC description to be useful 



So, are all KCMs (or more generally defect models) 
doomed to fail in describing thermodynamics? The point 
of view advocated in 0| is that after having appro- 
priately coarse-grained time and space a liquid can be 
strictly speaking considered as a KCM with trivial ther- 
modynamics and simple kinetic constraints that induce 
slow dynamics. As we have shown, this line of thought 
is clearly insufficient to recover thermodynamics. An- 
other point of view (see for example 0,1^) is to only 
assume that the effective dynamics for a lattice model of 
a glass-forming liquids is characterized by some kinetic 
constraints. In this case the trivial thermodynamics is 
not necessary but just a choice of simplicity. A quantita- 
tive model (arising from a real mapping from liquid dy- 
namics) would also contain effective interactions between 
particles and would therefore lead to non-trivial thermo- 
dynamics. One important fact remains though: even if 
the thermodynamics can be non trivial, it is only an indi- 
rect cause of slow dynamics, through the temperature de- 
pendence of the defect density, and not a driving cause of 



the slowing down. In our view, any eligible theory of the 
glassy state has to produce a convincing explanation of 
the remarkable connection between thermodynamics and 
dynamics, more precisely between the configurational en- 
tropy and the relaxation time |2^ . At the moment this 
seems to be an important missing piece in kcms. Hope- 
fully, future works on kcms and their possible general- 
isations will unveil whether or not some of these mod- 
els can pass this important test. A way to increase the 
entropy contribution of mobility defects and find a rela- 
tionship between dynamics and thermodynamics seems 
to be through extended defects. This is the path followed 
(using different arguments) in [Ullllllllllilliil (see 
also |30i. i31i| for quantitative computation of the excess 
configurational entropy in model systems). In all these 
scenarii, the presence of mobility regions, that are typi- 
cally interfaces (e.g. domain walls), is driven by a ther- 
modynamic mechanism with a (possibly avoided) critical 
point. 
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Addendum 

In the following reply. Chandler and Garrahan have 
amended their original model which now allows for two 
different species of excited cells - most of them are in this 
new version immobile, and very few are mobile. Even 
though the microscopic justification of this ad hoc as- 
sumption is unclear, the new fitting parameter, the ratio 
of excited mobile to excited immobile cells, indeed al- 
lows one to obtain both a large specific heat jump and 
a large relaxation time. In such a framework, however, 
the strong connection between thermodynamical and dy- 
namical fragility, which we emphasized above, becomes 
completely out of reach or very artificial (thermodynam- 
ics is dominated by immobile excited cells that are irrele- 
vant as far as dynamics is concerned). We also note that 
despite the introduction of another fitting parameter [the 
temperature dependence of the ratio of excited mobile to 
excited immobile cells] the theory proposed by Chandler 
and Garrahan still leads to a specific heat which has the 
wrong temperature dependence - note that the scale of 
the y-axis in their Fig. 1 makes it hard to see that their 
prediction for the excess of specific heat is in fact increas- 
ing with temperature, in contradiction with experiments 
that clearly show an opposite trend. It might be possi- 
ble, by introducing additional assumptions and fitting pa- 
rameters, to obtain yet another version of Chandler and 
Garrahan's model that "accounts for" these basic exper- 
imental facts. However, all this is at the expense of the 
simplicity, the predictive power and, more importantly, 
the very meaning of the underlying physical picture. 
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